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ABSTRACT 
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q \ Radio galaxies provide a means to determine the coordinate distance, the luminosity distance, the dimension- 
ed ■ less luminosity distance, or the angular size distance to sources with redshifts as large as two. Dimensionless 
coodinate distances for 55 supernovae and 20 radio galaxies are presented and discussed here. The ra- 
"si" \ dio galaxy results are consistent with those obtained using supernovae, suggesting that neither method is 
plagued by unknown systematic errors. The acceleration parameter q(z) and the expansion rate H(z) or 
j>. ■ dimensionless expansion rate E(z) can be determined directly from the data without having to make assump- 
t~-- tions regarding the nature or evolution of the "dark energy." The expansion rate E(z) can be determined 
| from the first derivative of the dimensionless coordinate distance, {dyjdz)~ x , and the acceleration parameter 
£Sj ■ can be determined from a combination of the first and second derivatives of the dimensionless coordinate 
,— 1 , distance. A model-independent determination of E(z) will allow the properties and redshift evolution of the 
| "dark energy" to be determined, and a model-independent determination of q(z) will allow the redshift at 
■ which the universe transitions from acceleration to deceleration to be determined directly. Determination 
q. of E(z) and q(z) may also elucidate possible systematic errors in the determinations of the dimensionless 
^ \ coordinate distances. 
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INTRODUCTION 

A primary goal of current cosmological studies is to determine whether the expansion of the universe is 
accelerating or decelerating at the present epoch, and what the acceleration or deceleration rate is. From 
C3 \ this, we can determine the global cosmological parameters, and study how structure evolved. 

One way to improve our understanding of the recent history of the universe, and begin to quantify the 
properties of the "Dark Energy," is to observationally determine the redshift at which the universe transitions 
from acceleration to deceleration if is is established that the expansion of the universe is accelerating at the 
present epoch. "Dark Energy" is discussed, for example, by Caldwell, Dave, & Steinhardt (1998), Turner & 
Riess (2002), Frieman, Huterer, Linder, & Turner (2002), Peebles & Ratra (2002), and Peebles (2002). 

Since astrophysics is a "spectator science," for which we can only collect data but cannot design and 
control experiments, unknown systematic errors are always a worry. The only way to convincingly establish 
a result is to have independent determinations that agree. 

THE BASICS 

The Robertson- Walker metric describes an expanding (or contracting) homogeneous and isotropic space- 
time, and has the line element 

dT 2 = dt 2 _ a 2 {t) ( + r 2 d0 2 + r 2 sin 2 Q (1) 
VI — KT J 

(see, for example, Weinberg 1972), where the cosmic scale factor a(t) is related to the source redshift z and 
the current value of the cosmic scale factor a Q : (a(t)/a ) = (1 + z)~ 1 . For light traveling from a source with 
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redshift z along a radial path, equation (1) implies that the coordinate distance to the source, a Q r can be 
obtained by integrating the equation 

dr/Vl - kr 2 = dt/a{t) . (2) 

In this universe, a source at redshift z, with intrinsic physical size D, and luminosity L, will have an observed 
angular size 9 given by 9 = D(l+z) / (a Q r) = D/cIa, where &a is the "angular size distance" d,A = (a r)/(l+z) 
and (a Q r) is the coordinate distance to the source. The flux / that could be detected from the source is 
given by / = L/[47r(a r) 2 (l + z) 2 } = L/[4ird1], where di is the "luminosity distance" di = (a r)(l + z) 
= (H )~ 1 y(z)(l + z); y(z) = H a r is the dimensionless coordinate distance (see, for example, Peebles 1993). 

Thus, if dA or dL to a source at redshift z can be observationally determined, then the coordinate distance 
(a a r) and the dimensionless coordinate distance y(z) to that redshift are known. 

The coordinate distance is related to the cosmological parameters through the equations a a J dr / y/l — kr 2 = 
f(a/d)dz, and (a/a) = H E(z). For k = 0, the coordinate distance is (a D r) = H^ 1 J dz/E(z). For a 
universe with components for which the equations of state Wi are time independent, such as a universe 
with quintessence (Caldwell, Dave, & Steinhardt 1998) E 2 (z) = ^^j(l + z) a \ where Wi = Pi/pi, and 
rii = 3(1 + Wi) for a component with non-evolving equation of state (see, for example, the Appendix of Daly 
& Guerra 2002). The deceleration parameter at the present epoch is q Q = —da /a = 0.5^2^(1 + 3wi), 
when Wi is time independent. 

Thus, one way to determine the cosmological parameters $7« and the equations of state Wi, is to determine 
the coordinate distance to high-redshift sources, then use equations given above, and solve for the cosmolog- 
ical parameters that yield the observed coordinate distances. If the equation of state is time-dependent, or 
if a rolling scalar field such as that proposed by Peebles & Ratra (1988) is considered, then these equations 
must be modified accordingly as discussed, for example, by Peebles & Ratra (2002). 

The cosmological parameters determined using this method then go into the equation for q a to determine 
whether the universe is accelerating or decelerating today. The equation for q(z), which is very similar to 
that for q a , then allows a determination of the redshift at which the acceleration is expected to go through 
zero, which marks the redshift at which the universe transitions from acceleration to deceleration. 

A method of using the observed coordinate distances to go directly to the acceleration parameter as a 
function of redshift will be described below. And, a method of using the data directly to determine the 
function E(z) without making any assumptions about the nature or redshift evolution of the Dark Energy 
will also be discussed. 

The beauty of using coordinate distance measurements to determine the cosmological parameters fij 
is that the true global cosmological parameters are determined; no corrections need to be made for the 
clustering properties of matter, so that contributions to Oj can not be missed or miscounted, and there are 
no biasing issues. 

THE RADIO GALAXY AND SUPERNOVA METHODS 

Two methods currently being used to constrain global cosmological parameters through the determination 
of the coordinate distance to high-redshift sources are the Type la Supernova method (e.g. Perlmutter et 
al. 1999; Riess et al. 1998), and the Type FRIIb Radio Galaxy method (e.g. Daly & Guerra 2002). 

For Type la Supernovae, there is one model parameter a, which goes into the determination of m e J^ , 
the effective apparent B band magnitude of the supernova at maximum brightness. This is related to the 
Hubble-constant free luminosity distance Dl via the equation 

m e J f (a) =M B + 5log(D L [ni, Wi }) . (3) 

If quintessence is being considered (e.g. Caldwell, Dave, & Steinhardt 1998), then Wi would represent the 
equation of state, and if an evolving scalar field is being considered, such as that proposed and studied by 
Peebles and Ratra (1988), then wi would be evolving with redshift (and would represent their parameter 

Mb is a constant obtained by fitting all of the data; it is related to the standard absolute magnitude of 
the peak brightness of a supernova Mg: M.® = Mb + 25 — 5 log(H ) (see Perlmutter et al. 1999). For k = 
0, and allowing for non-relativistic matter and quintessence, there are (N-4) degrees of freedom; this is also 
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the case for a universe with an evolving scalar field such as that discussed by Peebles and Ratra (1988), or 
a universe with space curvature, a cosmological constant, and non-relativistic matter. 

The Hubble-constant-free luminosity distance Dl is related to the dimensionless coordinate distance y(z) 
via the relation y{z) (1 + z) = Dl = H Q dL- 

For FRIIb radio galaxies, there is one model parameter f3, which goes into the determination of the ratio 
R* =< D > /D* (Daly 1994). This ratio also depends on the dimensionless luminosity distance Dl- 



k is a constant obtained by fitting all of the radio galaxy data. This fit also has (N-4) degrees of freedom. 
In these fits the dimensionless luminosity distance Dl, and the dimensionless coordinate distance y(z), is 
implicitly determined for each source, though it only factors out as a separate term when synchrotron cooling 
dominates over inverse Compton cooling with CMB photons in the radio bridge of the source, in which case 
R* =< D > /D* = (observables)(D L ) 9{l3) , where g{(3) = 3/7 + 2/3/3 (Guerra & Daly 1998). This is not a 
valid approximation for all of the sources in the sample. Here, this approximation has not been adopted, 
and an iterative technique has been used to determine the dimensionless coordinate distance to each source. 
Thus, these coordinate distances are valid for all of the sources in the sample. 

The dimensionless luminosity distance Dl and coordinate distance y(z) has been determined using equa- 
tion (3) for each source using the best fit value for Mb, obtained with the 54 supernovae in the "primary 
fit C" of Perlmutter et al. (1999) and the 1 high-redshift supernova published by Reiss et al. (2001), with 
the magnitude of this source corrected for gravitational lensing (Benitex et al. 2002). The dimensionless 
coordinate distance y(z) has been determined for each radio galaxy using equation (4) and the best fit values 
for k and (3 (along with their one sigma error bars) for the 20 radio galaxies presented by Guerra, Daly, & 
Wan (2000). These dimensionless coordinate distances are shown in Figures 1, 2, and 3. 

The dimensionless coordinate distances obtained for the radio galaxies are completely independent of 
those obtained for the supernovae. The conclusion from Figures 1, 2, and 3 is that there is good agreement 
between results obtained using Type la SN and FRIIb RG. This has also been demonstrated through detailed 
analyses. It is clearly the case when a cosmological constant and non-relativistic matter are considered 
(Guerra, Daly, and Wan 2000), when quintessence is considered (Daly & Guerra 2002), and when the 
evolving scalar field model of Peebles & Ratra (1988) is considered (Podariu, Daly, Mory, and Ratra 2002). 

The results shown here, obtained using 20 FRIIb radio galaxies obtained from the published literature 
and the VLA archive, will be improved and extended with 10 new sources that will be observed by O'Dea, 
Guerra, Daly, & Donahue at the VLA. 

Radio sources are observed out to very high redshift, and it would be easy to push this test to redshifts 
of three or four. The only thing that would be required would be a week or two of observing time at the 
VLA, followed by data analysis. This is not very demanding in terms of time, manpower, and funding. 

A MODEL-INDEPENDENT DETERMINATION OF q(z) 

The current method used to study the acceleration of the universe is to take the measured y(z) or Dl(z), 
determine best-fitting global cosmological parameters, and then use these global cosmological parameters 
to determine the acceleration of the universe as a function of redshift. In this process, assumptions must be 
made concerning the nature and redshift evolution of the mass-energy density of the "dark energy." 

A direct, empirical determination of the acceleration of the universe as a function of redshift can be 
determined using the data, without making any assumptions about the nature or evolution of the "dark 
energy." This can be done using the equation 



valid for k=0; if k 7^ 0, another term [kr(l + z)/(l — kr 2 )](dr/dz) must be added to the right hand side. 
Here, y is the dimensionless coordinate distance y = H (a r). 

Thus, equation (5) can be used to empirically determine the redshift at which the universe transitions 
from acceleration to deceleration without requiring assumptions regarding the nature and redshift evolution 
of the "dark energy." The supernova and radio galaxy data allow a determination of the dimensionless 
coordinate distance y to each source, at redshift z. These data can then be used to determine dy/dz, and 
d 2 y/dz 2 ; these can then be substituted into eq. (5) to determine q(z). 



R* (0,D L [n i ,w i ]) = K. 



(4) 



- q(z) = da/a 2 = 1 + (1 + z) {dy / dz)' 1 {d 2 y / dz 2 ) 



(5) 
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Fig. 1. Dimensionless coordinate distances y{z) to 20 radio galaxies and 55 supernovae as a function of (1+z). 
Note that the radio galaxy and supernovae determinations of y(z) are completely independent Radio galaxies 
are shown as open stars and supernovae are shown as solid circles. 
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Fig. 2. Focus on the low redshift end of Figure 1. Radio galaxies are shown as open stars and supernovae are 
shown as solid circles. 
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Fig. 3. The residuals between y(z) and those expected in a universe with Q m = 0.3 and Qa = 0.7, where y(z) 
is the dimensionless coordinate distance, shown as a function of (1+z). The error bars on the radio galaxies 
could be reduced by a factor of 1.4 as indicated on the figure if the source 3C427.1 is excluded and the reduced 
chi-squared is normalized to one, as discussed in detail by Podariu et al. 2003. Radio galaxies are shown as open 
stars and supernovae are shown as solid circles. 
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Equation (5) follows from the RW line element and the relation (1+z) = a /a(t). Our measurements of the 
coordinate distance (a Q r) move along the negative direction of dr, so eq. (2) with k=0 implies that a Q dr = 
— (1+z) dt, or (dz/dt) = — a~ 1 (l+z) (dr/dz)" 1 . Differentiating (1+z) = a /a(t) with respect to time implies 
that a = —a Q (l + z)~ 2 (dz/dt). Substituting in for (dz/dt), we find d = (1 + z) -1 (dr / dz) -1 . Differentiating 
again with respect to time, we find d = —(1 + z)~ 2 (dz/dt) (dr/dz)^ 1 [1 + (1 + z)(dr /dz) -1 (d 2 r /dz 2 )], 
which simplifies to eq. (5) using the expressions given here, and the relation y(z) = H (a r). Note, these 
expressions include the well-known equation (e.g. Peebles 1993; Weinberg 1972) 

H(z) = a/a = y/1 - kr 2 [d(a r) / dz}' 1 , (6) 
or, for k = and with H(z) = H Q E(z), 

E(z) = (dy/dz)- 1 . (7) 

Since eqs. (5), (6), and (7) are derived without any assumptions regarding the mass-energy components 
of the universe or their redshift evolution, they can be used to directly determine the function E(z), which 
contains important information on the "dark energy" and its redshift evolution, and to determine the 
dimensionless acceleration parameter q(z) directly from measurements of y(z). The use of the data to 
directly determine E(z) and q(z) is underway. 

DETAILS OF THE RADIO GALAXY METHOD 

The Radio Galaxy Method relies on a comparison of the average size of an FRIIb source determined by 
the mean size < D > of the full population at that redshift [< D >oc (a Q r)], and the average source size, D*, 
determined using a physical model that describes the evolution of the source. The two measures of the aver- 
age source size should track each other, so < D > /D* should remain constant, independent of redshift. The 
ratio depends upon observed quantities, the coordinate distance and the model parameter [3. When inverse 
Compton cooling with CMB photons is neglig ible, then =< D > /£>* = (Observables)(a r) 2l3 / 3+3 ^ 7 
(Guerra & Daly 1998). However, this is not a good approximation for all of the sources in the ra- 
dio galaxy sample. Therefore, this approximation has not been adopted here, and the full equation 
R* = k y^~ x ^l 7 (kiy~ 4 / 7 + fo)^/ 3-1 ), where k Q , k\, and £;2 represent directly observed quantities, has 
been used to solve for y(z) to each source. 

Details of the model are discussed by Daly & Guerra (2002), Guerra, Daly, & Wan (2000), Guerra & Daly 
(1998), and Daly (1994). The physics of FRIIb sources is discussed at length by Daly (2002). 

In brief, the average size of a given source is D* = f^t*, where i* is the total time that the AGN 
produces large-scale jets. Several different lines of argument, reviewed by Daly & Guerra (2002), suggest 

—3/3 

that oc L • , where Lj is the beam power of the source. If this relation is assumed, then oc 
(B L a L )- 2 ^ 3 v l ~ P/ * oc y-Wn+%1 7 ( fciy -4/7 + ^1-/3/3 _ For a typical value of p of L7; D ^ K (B L a L )- 1A v° L 4 . 

Note that the length of the source, D, does not enter into our determination of D*. Of course, D does 
enter into the determination of < D >. 

The determination of the most uncertain of these parameters, vl, can be studied using Chandra data, 
since Chandra data can be used to determine or constrain the ambient gas density n, and the radio galaxy 
model assumes that v\ oc B\/n. This is discussed below in the section "Tests of the Radio Galaxy Model." 

—3 /3 

The assumption that oc L- is consistent with models of jet production via the electromagnetic 
extraction of the spin energy from a rotating black hole. Daly Sz Guerra (2002) show that this rela- 
tion results if the magnetic field strength in the vicinity of the rotating black hole is given by B oc 
(a/mf/^-^ mW-^/^-W, or B oc (a/m) for (3 = 1.5, and B oc (a/m) 2 M a5 for (3 = 2, where a is 
the spin angular momentum per unit mass, m is the gravitational radius of the black hole, and M is the 
mass of the black hole (see Blandford 1990). 

The two other key assumptions of the radio galaxy model, that all FRIIb classical doubles at a given 
redshift have a similar maximum or average size so that the average size of a given source will be close to 
that of full population at that redshift, and that strong shock physics applies near the forward region of the 
radio source, have been tested empirically, and are consistent with the data. 
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TESTS OF THE RADIO GALAXY METHOD 

Three tests of the radio galaxy model and method are underway. These are tests using Chandra X-ray 
data; tests based on a detailed analysis of radio maps at multiple frequencies; and tests using a detailed 
comparison of radio maps with results from numerical simulations. 

Chandra X-ray data can be used to study the ambient gas density, which can be compared with that 
predicted. Donahue, Daly, and Horner (2002) have completed Chandra studies of 3C 280 (a radio galaxy) 
and 3C 254 (a radio loud quasar). The predicted ambient gas densities are consistent with the three sigma 
upper bounds placed by the Chandra data. Other radio predicted ambient gas densities will be compared 
with Chandra observations and bounds. 

Ten new RG with redshifts between 0.4 and 1.3 will be observed with the VLA in collaboration with Chris 
O'Dea, Eddie Guerra, and Megan Donahue. These high-resolution observations will allow several tests of 
the RG model, and will be compared with detailed numerical simulations tailored to match the properties 
of these sources, which is being done in collaboration with Joel Calvalho and Chris O'Dea. These studies 
will be based on the numerical work of Calvalho & O'Dea (2002a, b), and will help to identify the physical 
processes that must be accounted for the in modeling FRIIb sources. 

In addition, complete radio galaxy samples of FRIIb sources that go to redshifts of three or four will be 
investigated. These could serve as the parent populations in the application of the RG method to redshifts 
greater than two. This could take us to very high redshift very quickly since radio observations are relatively 
quick, easy, and do not require major new instruments. 

CONCLUSIONS 

In a spatially flat universe with non-relativistic matter and quintessence, radio galaxies alone indicate 
with 84% confidence that the universe is accelerating in its expansion at the present epoch (Daly & Guerra 
2002). Results obtained using the Radio Galaxy Method out to redshifts of two are consistent with those 
obtained using the Supernova Method out to redshifts of about one. 

A model-independent way to use the supernova and radio galaxy data to determine the dimensionless 
expansion rate E{z) and acceleration parameter q(z) is presented and discussed. A direct determination of 
q(z) that is independent of assumptions regarding the nature and evolution of the "dark energy" would allow 
the transition redshift from acceleration to deceleration to be determined, and would help to identify any 
systematic errors that might plague either the radio galaxy or supernova methods. A direct determination 
of E(z) would help to quantify the properties and redshift evolution of the "dark energy" and to identify 
potential systematic errors in the methods used to determine y(z). 
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